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Abstract 



The Higgs couplings to matter fields are proportional to their masses. Thus 
Higgs amplitudes can be obtained by differentiating amplitudes without 
Higgs with respect to masses. We show how this well-known statement can be 
extended to higher order when renormalization effects are taken into account. 
We establish the connection with the Callan-Symanzik and renormalization 
group equations and consider also pseudoscalar Higgs couplings to fermions. 
Furthermore, we address the case where the Higgs couples to a heavy particle 
that is integrated out from the low-energy effective Lagrangian. We derive ef- 
fective interactions where mass logarithms are resummed by renormalization- 
group methods, and give expansions of the results up to next-to-leading order. 



1 Introduction 



The Higgs boson provides a simple mechanism to accommodate massive vec- 
tor bosons and fermions in the standard model. Present-day experiments, 
together with calculations of Higgs interactions up to two-loop order in some 
cases, have been used to establish lower limits for Higgs masses approaching 
the W and Z mass range, and the next generation of colliders may give a 
definite answer to the question of its existence. 

A widely used tool in the study of Higgs interactions are low-energy the- 
orems (soft-Higgs theorems), which play a role comparable to the low-energy 
theorems for pion amplitudes in hadronic physics. They rely on the fact 
that the explicit breaking of scale invariance by the Higgs interactions can 
be employed to relate tree-level amplitudes with different numbers of zero- 
momentum Higgs fields [|IJ. This theorem has been extended to one-loop 
amplitudes, where scale invariance becomes anomalously broken, and it has 
been observed that there exists some connection with the scaling functions 
(beta functions) of renormalization group theory ||. Various applications 
can be found in the literature H M, and recently with its help two- loop am- 
plitudes were calculated in the heavy-top limit 0, |], [7], || , where algorithms 
were devised to take into account the renormalization effects. However, the 
role of scale anomalies in higher order has remained unclear in the present 
context H , and thus the precise form of the theorem for renormalized ampli- 
tudes in the general case has remained unknown. The purpose of the present 
paper is to clarify this issue and to provide a general survey of soft-Higgs 
theorems in the framework of renormalized perturbation theory. 

Since it is conceptually simpler, we shall discuss first the soft-Higgs the- 
orem for pseudoscalar Higgs bosons, which exist in models with an extended 



Higgs sector. In that case the relevant symmetry is chiral invariance |10 
anomalously broken by the well-known triangle anomaly JTI]. Next, the the- 
orem for scalar bosons will be developed, where the anomalies proliferate. 
We shall show how they are controlled by the Callan-Symanzik equation, 
and give the explicit form of the theorem both in on-shell and minimal sub- 
traction (MS or MS) schemes. The latter allows for the introduction of 
effective-theory methods, which already have been applied in ||. That the 
effective-theory picture is appropriate, follows from the observation that the 
Higgs coupling to the heaviest particle (e.g., the top quark) is dominant, 
and since the soft-Higgs theorem applies at low energies (and small Higgs 
masses), such a particle should be integrated out from the low-energy the- 
ory. Furthermore, when this method is used, logarithms of large mass ratios 
are easily summed by the renormalization group. We shall derive the form of 
the soft-Higgs theorem in the effective theory and show that this framework 
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provides a natural description of all coefficients in terms of scaling functions, 
which we shall give in some detail. In an appendix we give the formulas in a 
form which is directly applicable in a next-to-leading order calculation, and 
demonstrate their use in a sample calculation that can be compared to the 
calculational methods used in the literature. 



2 Pseudoscalars and chiral symmetry 

Before we consider the scalar Higgs, let us investigate the couplings of a 
pseudoscalar (CP-odd) Higgs boson, predicted, e.g., by supersymmetric ex- 
tensions of the standard model, to fermions. To keep things simple, we allow 
only one external pseudoscalar in the amplitude. 

We apply an infinitesimal global chiral transformation 

Sip = -ieg A j 5 ip/2, (1) 

where g A is the coupling of the fermion ip to the pseudoscalar Higgs A: 

£a = gA&lsMi ( 2 ) 

The variation of the fermion mass term introduces an additional term in the 
Lagrangian 

5C = eg A m^i^ 5 ij;, (3) 
which may be interpreted as an imaginary contribution to the masses of left- 



and right-handed fermions [10|. To first order in e, the fermion propagator 
is modified into 

n ( , v .^ + m + ieg A mj 5 m 

G e {k) = i . (4) 

k A — mr 

Thus we have the relation 

¥e C = ™— = SC. (5) 

In addition to @, the variation (HD introduces terms mixing the scalar (H) 
and pseudoscalar (A) couplings to fermions, since these interactions also 
break chiral invariance. This mixing can be compensated by a redefinition of 
the H and A fields (in the light-Higgs limit) and will not be considered here. 

Apart from its coupling to fermions, the A field can couple to other Higgs 
or Goldstone fields. The corresponding interaction operator will be called A. 
This term is not obtained by the differentiation with respect to e: 

5C = m^- - A. (6) 

5A v 1 
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For the full amplitude, differentiation with respect to e obviously will give 
no diagrams with A insertions. These have to be calculated separately, but 
in the limit of light A they can be neglected, since the Higgs masses are 
proportional to self-couplings. 

The relation (U) may be extended to a relation for T, the generating func- 
tional of one-particle irreducible (1PI) vertices, which is in tree approximation 
equal to the action / L. First we write down the identity describing the ir- 
reducible interactions of one zero-momentum A in terms of zero-momentum 
operator insertions: 



I Ja T{0)1h=a=0 = I 5C ^ r(0) l^=o + / * 4 r (0) \h=a=o. (7) 



(We use A Jj T as a shorthand for the 1PI generating functional with an 
insertion of the renormalized operator A.) It holds in the presence of quan- 
tum effects, up to scheme-dependent universal corrections (denoted by a a) 
referring to the particular renormalization conditions imposed on the A field 
(cf. Sec. |): 

m(l + a A ) J jjT\ H=A= v = j ' 8C^ V\ H=A=a + j A$ T\ H=A=0 . (8) 

This is seen by investigating the diagrams contributing to both sides. 

The nontrivial part is the quantum extension of the Ward identity (EJ) of 



chiral symmetry. However, the answer is well known ||11|| : 

^r(e)| e=0 = JSC Y - sg A fix [fP] ^ T. (9) 

where s is a constant equal to its one-loop value, and F (F) denotes the 
(dual) field strength tensor. Thus we obtain 

mil + a A ) J — r = — r + sg A J tr [fP] ^T + Ja^T, (10) 

where e, H, and A have to be set to zero after differentiation. 

The relation ( |I0~D is readily verified in an explicit calculation. If the 
integrand contains an open fermion line, the insertion of a zero-momentum 
A gives the expression 

■ ft + m tf + m . ig A ls Id 

l T2 2 l ^ Al ^) l T2 2 = % T2 2 = —TrGe(k)\e=0, (H) 

k — m A k z — m z k A — m z m oe 

where we have used the anticommuting nature of 75. Thus the theorem is 
trivially satisfied. In fermion loops, the behavior of 75 is accounted for by 
the triangle anomaly term in flT0|). 
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Since we work only to first order in e, it appears only in the numerator 
and the differentiation does not affect the large-momentum behavior. Thus 
there is no room for renormalization to introduce further anomalies. If we 
were to derive the interaction of several pseudoscalar particles, apart from 
other complications the quadratic term in e would give a mass shift, so that 
renormalization effects came into play. Then the scale anomalies introduced 
in the next section would have also to be considered. They are CP-even and 
therefore do not occur in the first-order A couplings. 

As mentioned in the introduction, the theorem ([II]) is valid for vanishing 
four-momentum p&. Only when the mass can be neglected compared to 
m — in which case also the contribution of A will be small — this approaches 
the amplitude for an on-shell A particle with p& — 0. For this reason, (ID) 



may be called a soft-v4 theorem (for applications, see ref. ||, and references 
cited therein). 

3 Scalars and scale transformations 

Turning over to scalar Higgs couplings in a theory such as the Standard 
Model, we divide the Lagrangian into three parts 

£ = C + C GF + C H , (12) 

where the first term denotes the Lagrangian of gauge boson and matter 
fields, including interactions with the Higgs fields, the second part is the 
gauge-fixing and ghost part of the Lagrangian, and the last term contains 
the pure Higgs Lagrangian, i.e., the self- interactions of the physical Higgs H 
and the Goldstone fields. 

We assume that mass terms in £ are generated only through the Higgs 
vacuum expectation value v, so that 

£ {m, H) = Coirrii + gji) (13) 

where the mass is given by 

rrii = giv. (14) 

Similar to the connection between pseudoscalar couplings and chiral transfor- 
mations described in the preceding section, scalar Higgs couplings are related 
to scale transformations. These act on the fields as 

Sip = (d^p + x ■ d)tjj, (15) 
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where is the canonical dimension of the generic field ip. With = 3/2 
for fermions and d^ = 1 for bosons, the variation of the action reads 

SC = I (v^ + 5C GF + 8CrA , (16) 

since scale invariance is broken by a nonvanishing value of v, and by dimen- 
sioned parameters in the gauge-fixing and Higgs parts of the Lagrangian. 



Using (13), this can be rewritten as 




or 



M = I ( V TTr - V J^ + ^GF + SC H ^j (17) 

Similar to the operator A in the preceding section, the operator TL sum- 
marizes Higgs self-couplings and couplings to gauge fields generated by the 
gauge fixing. The diagrams with TL insertions have to be calculated explic- 
itly. When they contain Higgs (Goldstone) self-couplings, in the light-Higgs 
limit they can be neglected compared to diagrams with couplings to heavy 
particles. However, the mass parameter in the gauge-fixing part (in a 'tHooft 
gauge, for instance) introduces a spurious variation 5Cgf that is not related 
to a Higgs coupling. In one-loop order it can be separated just by omitting 
the derivative with respect to this gauge-fixing mass in the relations below, 
but from the two-loop order on it becomes entangled into the renormaliza- 
tion procedure. There are several possibilities to deal with this complication: 
one could either use the background-field method employed in fl2|| , which 



requires the calculation of more diagrams, or impose the additional restric- 
tion that also the gauge boson masses have to be neglected, or turn over to 
a mass-independent renormalization scheme. The latter will be done in the 
next section. 

The extension of (|17T) to irreducible vertices has the form 



(1 + a H )v J -^T = J 5C\^T + ju^T. (19) 

The unrenormalized diagrams which do not involve Higgs (Goldstone) self- 
couplings, contributing to the relation (|19|) , are identical on both sides. The 
same is true for the counterterms, except for the renormalization of the exter- 
nal Higgs field: In the on-shell scheme the left-hand side is renormalized with 
the momentum pn flowing into the vertex satisfying p 2 H = m 2 H) whereas the 
right-hand side is defined for vanishing momentum pn- Corrections propor- 
tial to m 2 H are consistently absorbed into TC, so that there remains an overall 
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multiplicative factor, denoted by 1 + an- Any possible correction involving 
Higgs (Goldstone) self-couplings is also absorbed in 7i. 

We can employ the Ward identity of scale invariance to get rid of the 
operator 5C in ([19]) and obtain in this way a nontrivial statement. Although 
scale invariance is broken by quantum corrections ||13||, the quantum action 



principles [|14[] tell us that all corrections (anomalies) are given by a linear 
combination of local operator insertions with the appropriate quantum num- 
bers and dimension. The resulting anomalous Ward identity of scale trans- 
formations is known as the Callan-Symanzik equation For an abelian 



Higgs model a detailed derivation is given in Here we simply generalize 
the result for a general Higgs model in an on-shell scheme: The relation (|19|) 
can be replaced by 

+ J Hl^T^'-l (20) 
for an irreducible vertex with external fields of species i, which is given by 

r(»i-) — r (21) 

The £i are gauge parameters, and Tt can be neglected, as discussed above. 

The coefficients a, (3, 7, 5 have to be determined order by order in per- 
turbation theory. They are universal for all processes, but they are mass- 
dependent and are not simply related to the familiar scaling functions in 
a mass-independent scheme, which will be introduced in the next section. 
Some of them are determined by the symmetry. In particular, all but one 
of the coupling constants are usually expressed in terms of masses, so that 



their beta functions may be identified as mass beta functions |12| , due to the 
normalization conditions one imposes. 

Introducing the concept of bare parameters, which can also be used for 



a simple derivation of the Callan-Symanzik equation |13| , the coefficients 
can be expressed as derivatives of the renormalized parameters with respect 
to bare masses. The relation (^) thus involves the operations carried out 
in [[], [7J in reverse order: If instead of renormalizing after taking the derivative 
with respect to bare masses, the renormalized amplitudes are differentiated, 
one has to correct for the derivatives of the renormalized parameters. These 



are the corrections summarized in fl2"0"|) . 
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4 Renormalization group 



The relation fl20|) has been derived within the context of an on-shell scheme, 
where all dimensionful parameters are expressed in terms of the physical 
masses of the theory. When a minimal subtraction scheme (MS or MS) is 
used, the derivation is greatly simplified. The price are complicated, but 
calculable, relations of the parameters to physical observables. 

Taken literally, a minimal subtraction scheme does not subtract the tad- 
pole diagrams completely. This is a minor problem: In terms of irreducible 
diagrams, only the Higgs one-point function is affected, and since it is a con- 
stant, the finite part can completely be subtracted without affecting prop- 
erties of the dimensional renormalization procedure that follow from mass 
independence: tadpoles are simply omitted. 

To lowest order (tree level), the relation ( |20| ) reads 

v [ — r< ni --) = ^m.-^-rt™ 1 '-) + In \ r( ni >~\ (22) 

J 5H OTflj J 

Here we exclude the mass term in the gauge fixing from the derivative, so 
that 7i does not contain spurious couplings and is truly negligible in the 
light-Higgs limit. 

In the dimensional renormalization scheme, the quantum action principles 



hold in the strong sense [16], so that (22) is valid without corrections even 



after renormalization, if the running parameters at the scale fi are inserted 
everywhere. Stated differently, the coefficients [3 and 7 appearing in (p0|), 
which are related to mass derivatives of the renormalized parameters, vanish 
identically: the renormalization factors are mass-independent. 

It is now easy to get rid of the mass derivatives in favor of familiar renor- 
malization group coefficients, if we have a one-scale problem. The case of 
multiple scales will be considered in the next sction. We thus put all other 
masses and external momenta to zero, ignoring infrared divergences. (If some 
care in the renormalization of sub divergences is taken, they may be regulated 
dimensionally.) We use the MS or MS scheme. Then any renormalized vertex 
function with mass dimension d and rig external light fields of species I can 
be expressed as 

r(m(/i), m) = m(fx) d f (In g(ji)), (23) 
with a dimensionless function T. Defining the mass anomalous dimension 

d\nm(u) , . , . 

» — = -7m(A* , 24 
am// 
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we derive 



d -I» = (d-(l + lm ) - 9 . - + ftA)r(mfc),/*)- (25) 



din// \ (91nm(yu) 

On the other hand, the vertex function satisfies a renormalization group 
equation 

ik rw = -? 7 ' r(rt ' (26) 

so that 

° r(m(n), n) = — L_ f A A + d + !?0 r(m(/x), //). (27) 



dlnm(fi) ' 1 + 7 m \ % 2 

The definitions of the scaling coefficients are 

A = ~7T~ 9ii (28) 

dm// 

% = 7T- lnZ - ( 29 ) 
am/i 

Thus the soft-Higgs theorem in the dimensional renormalization scheme has 
the form 



v 



i_r(m,-) 
5H 



1 + 7m V 9 9i 2 J 

+ Jn\^ r( fti ->. (30) 



This looks very similar to (|20f) , but without the derivative with respect to the 
mass: In this scheme the coupling of the Higgs is given by scaling coefficients 
only, for vanishing external momenta. 



5 Heavy and light fields: effective theory 



The arguments leading to (p0|) are sufficient in a single-scale problem, when 
the scale \x can be chosen not much different from m, and infrared divergences 
can be ignored. However, in particular when QCD corrections are consid- 
ered, it is customary to resum logarithms because of the comparatively large 
coupling constant. In this situation it is mandatory to apply the method of 
effective field theory |[L7|| , since below a mass threshold the renormalization 
group of a mass-independent scheme is not able to resum logarithms correctly 
(see [[TH] for a discussion of this point). 
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This framework is also helpful to understand the role of different scales 
in the presence of light masses and small momenta. It clearly separates 
the infrared behaviour from the ultraviolet, and all coefficients in the soft- 
Higgs theorem can be expressed in terms of anomalous dimensions and beta 
functions. 

The heavy particle is integrated out at a scale //o of the order of its mass 
rrih, with the effect that the Lagrangian £ which can be divided into a part 
£i that contains only light fields, and the rest £h 

£ = £ e + £ h , (31) 

is replaced by an effective Lagrangian 

C = £ e + Y,CiOi. (32) 

i 

The operators Oi are of increasing dimension, divided by appropriate powers 
of rrth- They consist only of light fields. In particular, they contain operators 
of dimension four or less which are already present in It is convenient to 
absorb those into a redefinition of parameters: 

£ = £ e + Yl £A- (33) 

dim>4 

The ordering of operators according to their dimension is appropriate when 
no light masses are around. Otherwise it simplifies the discussion if the series 
is organized in terms of powers of instead. 

One should keep in mind that the limit under consideration is not exactly 
a heavy- mass limit for, e.g., the top quark (m t — > oo would imply strong 
coupling), but rather a small- coupling limit for the Higgs self-coupling. The 
use of effective field theory methods in this limit for the resummation of 
logarithms is somewhat unusual, but in perturbation theory where the mass 
of a particle and its coupling are clearly separated, there is no real difference 
to an ordinary theory with large mass ratios. In particular, the fact that 
a particle like the top quark may be of a non-decoupling nature introduces 
no practical difficulties, since we stay in the weak-coupling regime where the 
mass is still small compared to the strong-coupling scale Aw. 

The coefficients Ci are determined as power series in the couplings by a 
matching calculation. It involves calculating the difference of diagrams in 
the effective and full theories, where the mass has to be kept in the full- 
theory diagrams. With respect to the heavy particle it thus incorporates 
the change from a mass-independent to essentially an on-shell scheme (for 
details, see e.g. |1|)- This is necessary to ensure that below threshold the 
correct logarithms are summed. 
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An important property of the matching coefficients Cj in ( p2| ) is that they 
are infrared safe quantities: They are infrared convergent, and they contain 
no logarithms of light parameters me. Thus any dimensionless matching 
contribution A is a function of the ratio of the two dimensionful parameters 
m h{Ho) an d Ho, which appear only logarithmically, and of the dimensionless 
couplings g(/J,o): 

A = A ( ln ^fc)^M). (34) 

It is independent of light masses and momenta, up to power corrections 
0{me/mh,pf/mh) which are absorbed in the higher-order terms. 

Since masses are renormalized multiplicatively, the vector d/dm(fi) scales 
contravariant to the mass vector m(/x), and the scalar product which appears 
in fl22|) is invariant: 

? m ' M 9^) = ? m ' ( "° ) ssk)- (35) 

Thus we can take derivatives at the matching scale where the heavy particle 
is removed from the theory. In the effective theory, the heavy mass mh(no) 
appears implicitly in the coefficients, whereas the light masses rrie(fio) (reex- 
pressed in terms of effective masses rhe(n)) remain dynamical parameters. 



6 Mass dependence of effective parameters 

Before we state the soft-Higgs theorem in the effective theory, we first have to 
calculate the dependence of the parameters in the low-energy effective theory 
on the heavy mass m^. Let us consider first the running coupling constants. 
In the full theory, they satisfy 

</(//) = /%(//)) (36) 



cHn/1 

where g and (3 are vectors so that this is a coupled system of differential 
equations in the general case. By definition, in a mass-independent scheme 
such as MS the value of a running coupling constant does not depend on 
masses 

Tj~r~~ gin) = 0, (37) 
a In m/j 

although the existence of the heavy particle affects the beta function. 
The transition from fl32|) to fl33|) is reflected in the matching condition 



go = gM = gM + A^in jl, ^M), (38) 
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where Ag is a polynomial function of g starting with the cubic term. Only 
diagrams containing the heavy particle contribute to Ag. In the effective 
theory, the running coupling constants satisfy 

d 



dliafi 



(39) 



where (3 is obtained from j3 by omitting the diagrams containing the heavy 
particle. The renormalization group can be used to resum logarithms and to 
obtain the value of g at another scale \x 



g{n) = g(\n-^,g ). 

When evaluated exactly, g(fi) is in fact independent of fi : 

dg(^) = Q 
d In /i 

Inserting the definitions (40) and (|38|) into this identity, and using (|3( 
we derive the relation 

dg(p) 1 



(40) 

(41) 
and 



d\nm h (no) 



1 + 1hh(fM>) + 



dg 



1 + 



dAgjfMp 
dg(^o) 



(42) 



The other parameters of the effective theory include the matrix of field 
renormalization factors ( (with Z = ( T ( being the coefficient of the kinetic 
term in the effective Lagrangian), the masses rhi of light fields (we consider 
only fermions), and the coefficients Ck of additional operators of the order 
In the full theory, the renormalization group equations are 



GHn/i 

cHn/i 
dnih 



1 



C0)70), 



c?ln/i 

The matching conditions yield 



-m h {n) 7w(A*) - ™>i(iJ,) 7«(//), 



CM 



CM [l + AC(ln 
m h (Ho) Aa(ln 



Mo 



> 9(1*0)) 



mh(Mo) 



1 + Ar(ln 
AC fe (ln-^ 



Mo 



**- <?M) 



#M) 



(43) 
(44) 
(45) 

(46) 

(47) 
(48) 
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where ( is equal to ( projected onto the space of light fields. 

In the effective theory, the coefficients satisfy renormalization group equa- 
tions 



drhe 
dlnfi 



2^)7^), 
-m e (fi)^f ee (fi), 



and 



d 



+ ll C k {fi) + £ iJ lh ...CM CeM • • • = 0. 



(49) 
(50) 

(51) 



\d\nfi , T,£=k 

By definition, operators mix among each other only if they have the same 
power m^ k as prefactor. The nonlinear terms appear for k > 1; they origi- 
nate from time-ordered products of lower- dimensional operators mixing into 
the local operators. The renormalization group equations can be solved iter- 
atively for increasing dimension so that the nonlinear terms play the role of 
the driving term in an inhomogeneous differential equation for the coefficients 
with index k. 

The solution of the renormalization group equations can be cast into the 
form 



C(fi) = C(fi )R(lnf o ,g(fi )), 
m e (fi) 

cm 



■m e (fi ) R u (}n^,g(fi )) 



R k{^f o ,g(fi Q ))C k (fi Z 



+ £ Rj lh ...(ln^,g(fi ))C e MC e2 (fi ) 

Y>l=k 

Requiring fi independence of these expressions, we obtain 



(52) 
(53) 

(54) 



c 1 



d In rrih 



1 + 7w» 



L_{«- 1(1+ Acr^(i + Ac)«-^ 







+ R-t(l + A()- 1 (3— (1 + AQR 



dg 



(55) 



drhi 
d\nm h 



d 



AaR u + p— (AaRu 
dg V 



- 7w(l + At)R u + Act Rulu 
d 



P-((l + Ar)R u ) 
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7«(1 + At)R u + (1 + At)R u % 



(56) 



dC, 



d 



(57) 



In the first equation, 7 is the anomalous dimension matrix 7 projected onto 
the space of light fields. All full-theory and matching coefficients are evalu- 
ated at the scale /z , and the effective theory coefficients (denoted by a hat) 
at the scale \i. It is convenient to choose 



Ho = m h (fjto) = m h (m h ), 



(58) 



so that all logarithms vanish in the matching coefficients. 

When in the effective theory the logarithms have been resummed, the 
solutions R x (n) are usually available as functions 



R x (/i) = R x (g ,g(n)), 



(59) 



with g from ([38|), containing no explicit logarithms. Then the derivative 
with respect to g(/U ) is given by 



dR x 
dg(n c 



PM [1 + 
I /3(/io)fl4 



dAg \ dR x 



dg(no)J dg 

dAg \ dg(n) 
d9(j*Q) ) d 9o 



(60) 



dR x 



7 Soft-Higgs theorem in the effective theory 

We are now in a position to derive the soft-Higgs theorem for the effective 
theory. Applying the chain rule to (|22|), we obtain 



d 



t mi(l+Ar) L 



+ 



+ 



drhf d 



+ 



dg d 



d In m/j drhi d In rrih 
' dC k 



-1 dQ 



dining 



din rrih 



o k + n) ^r< ni -\ 



p(m,...) 



(61) 
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with Q = [Tii) 1 ! 2 being the normalization of the light particle i which appears 
ri£ times in the vertex function. (We have neglected the mixing in the light 
sector for simplicity.) It is understood that Yukawa couplings are expressed 
in terms of masses only after the differentiations have been carried out. 

This equation replaces fl30f ) in the general case, taking now full account 
of light masses and nonvanishing momenta. When the explicit expressions 
(|42|), (|55|), fl56|) , and ( |57|) are inserted, a simple pattern emerges: The coef- 
ficients essentially consist of differences of scaling functions (beta functions 
and anomalous dimensions) in the full and effective theories. In addition, 
beginning from second order the mass dependence hidden in the matching 
contributions has to be accounted for. (In the non-decoupling case, where 
heavy and light masses are allowed to mix, a matching contribution enters 
already at leading order.) This generalizes the observation in || that the 
leading contribution to, e.g., the H — ► 77 decay amplitude, is determined 
by the contribution of heavy particles to the beta function, which is equal to 
the difference of beta functions in the full and effective theories. 

The parameter set of the effective theory (denoted by a hat) is reduced. 
Along with the heavy mass it does no longer contain the couplings of the 
heavy particle. Thus there are no problems arising from the fact that its 
Higgs coupling is related to its mass and should be defined at the same 
scale. If we were to use the full renormalization group in a mass-dependent 
scheme (to account for mass effects), with a dynamical heavy particle below 
threshold, we would run into difficulties EO . 



8 Conclusions 

The soft-Higgs theorem has found a wide range of applications in standard 
model calculations, for Higgs amplitudes in the limit of small Higgs mass 
and momentum, where it can be used for approximations and as a nontrivial 
check for complete calculations. The aim of this paper was to clarify its 
meaning in the context of renormalized perturbation theory. We have shown 
that there is a close connection to the Ward identity of scale invariance, the 
Callan-Symanzik equation. Similarly, broken chiral symmetry is related to 
the coupling of pseudoscalar Higgs fields. 

The exact form of the renormalized soft-Higgs theorem depends on the 
renormalization scheme one has imposed. We presented results in the on-shell 
scheme, and in a minimal subtraction scheme (MS or MS), where things sim- 
plify, and in particular the coefficients are mass independent. If large mass 
ratios are present, so that logarithms have to be summed up, heavy particles 
have to be integrated out below threshold, and the soft-Higgs theorem as- 
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sumes a different form. We have calculated the coefficients in a fairly general 
way, so that it should be straightforward to use them in a particular problem. 
The many possible extensions of the standard model leave plenty of room for 
new particles and interactions up to the TeV range, where the soft-Higgs 
theorem could be valuable in the calculation of Higgs interactions. 
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Appendix: NLO expansion 

To establish the connection to applications of the soft-Higgs theorem that 
have been considered in the literature, we expand the various terms up to 
next-to-leading order in the coupling constants. The complicated structure of 
the full standard model, including QCD, forces us to maintain full generality 
and to keep the mixing of the various coupling constants, fields, and masses. 
However, there may be no need to resum logarithms, as it has been the case 
in existing applications (for instance, with the present knowledge of Higgs 
and top masses, ln(mt/m#) is not particularly large). For this reason, we 
neglect higher-order logarithmic terms in the formulas below and express 
everything in terms of g % = g l (^o), ( = C(^o), and me = me(fio), where 
the matching point fi is chosen equal to m/ l (/i )- With = g l gi ■ ■ -, and 
a summation convention for upper indices understood (all tensors may be 
chosen symmetrically), we expand up to next-to-leading order: 



CO*) 

rhje(ii) 
CM 



j + g ^(ffilnJL + A gg), (62) 
V m h J 

C(H //" A,;') /,'(//). (63) 

m t (l + g ij Ari j ) + m h (g'-'Acr- 1 + g^Aaf 1 )] R u {n), (64) 



ml 



1 " 9 ij lli ^ 



g lm "-ACl 



Ira- 
ki 



(65) 
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R(/i) and Ru(fJ>) are chosen to be solutions of the renormalization group 
equation valid to order g A . The initial Wilson coefficients AC& are assumed 
to be of order g n , with n > 2. The expansions of the scaling functions are 



tig) = g ijk tii k + g ijklm tii2 lm , 
i{g) = g ij ^ + g ijkl if\ 



and so on. A straightforward calculation then gives 



dg l 



dlnrrih 



r 1 



dlnrrih 



drhe 
dlnrrih 



dC, 



din rrih 



ijk 



g l 



ijklm ( j^3 klm 



i2 



njklm\ 
Pi2 ' 



J 



(66) 
(67) 



(68) 



+ g tJklm (Pii - Mi) 

_|_ 2g i 3 klm 



Im 



l,lj h 



Ihhl 



tia ~ tin) In ^ + AgSfi? - 0>}Ag$ 



1 



u'-A 
+ g 



2 

ijkl 



ii j ) 



(69) 



I( 7 f'-^)-I 7 ^( 7 f-^) 



+ 7? (Ai - til) In ^ + 2A^ti! - li^gli 
z 4 m>h 



ki 



= m h {g ij {^- 1 Z 1 ) 



(70) 



ijkl 



ijkl 
lU2 



+ g ijkl [(Act. 

- (A(7? - 7^1) 7m In ^ - 7 &i (A^f - 7S1) 
+ A^ 7m - 7S1A7* + 2AaM] } 
+ m< I g lJ ( 7m - 7m ) 



+ g 



ijkl 



{lm lm ) lihi {lui 7m) 



- (7& - 7m) 7m In f - - -& (Aa* - T&) 



+ A^ 7m - 7m Arf + 2At?(% 
+ 27m (#? - ffi) ln-^ + 2 7m A^ 

_L {_ fcff im- _ ( 7 ^ + n^flg'} C£T 
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m 



(71) 



Here the bracket [•, •] denotes the commutator of two matrices. 

These relations are complicated. In practical applications, one usually 
looks for higher-order effects only in the QCD coupling, so that many terms 
can be dropped. Also the resummation of logarithms is then simplified. 



Example 

To demonstrate the method, we take the process H — ► bb which has been 
investigated in ref. ||. Let us consider the contribution of the diagram 
in Fig. [I] in the MS scheme, together with its counterterms (Fig. ^j. By 
A (B,...E) we denote the residue of the 1/e pole, where the space-time 
dimension is D = 4 — e. We consider only the vectorial part of these self- 
energy diagrams, which is proportional to fb- The analysis of the parts 
proportional to m&, and of the other diagrams contributing in the given 
order, proceeds along similar lines. 

The relevant formula is Three coupling constants are around: g^, gt 
(Yukawa couplings), and g s (QCD coupling). However, terms proportional to 
gb do not contribute to the wave-function renormalization in leading order, so 
that all expressions in our example are proportional to g 2 g 2 . The counterterm 
diagrams correspond to renormalizations of g s and g t . Of course, the net 
renormalization of g s is zero here, due to the QCD Ward identities, but the 
compensating terms are provided by different diagrams. 

First, we collect the terms which are themselves of order g 2 g 2 and have 



to be inserted into (^H]) at tree level. The following terms are nonvanishing: 



g 2 s gh* ss /2 = (2A-2B-2C)/2, (72) 
-<^ f 2 7f72 = ~(2D)/2, (73) 
^ 2 ACfAT = C/2. (74) 

For the calculation of 72 we had to subtract the counterterms for the UV 
divergent subdiagrams (B and C) twice. The effective theory term 72 con- 
tributes one diagram of order g 2 , with an insertion of the matching contribu- 
tion proportional to g 2 . The matching term A£ is proportional to g 2 and has 
to be multiplied with one contribution to the g t beta function; this is easily 
seen to be again equal to C/2. The prefactor in the left-hand side of ([74] ) 
is 1, not 2, since 2A£i/3i has to be distributed among Fig. [I] and its mirror 
image. 

Furthermore, we have to consider a contribution of ( |68|) in order g 2 , in- 
serted into the master equation (pTl) and evaluated in the effective theory to 
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order g 2 s . This is equal to 



AfW = 5. (75) 

Taking everything together, the sum of the terms pertaining to Fig. [I], in- 
serted into (^) with n& = 2, is 

2A- B -C -2D, (76) 

We have assumed that all masses are kept nonzero. 

This may be compared with the approach of ref. 0. There the light 
masses are neglegted, and the resulting IR divergences are regulated dimen- 
sionally. This implies that B and D are identically zero, and their contribu- 
tions are absorbed into A. The calculation of Fig. |l|, and the subsequent dif- 
ferentiation with respect to m" gives a contribution 2A. The renormalization 
of g s and g t (the former being in fact obsolete) is equivalent to subtracting 
B and C once. Thus we have again the result 

2A- B -C -2D. (77) 

If the on-shell scheme is used, E is also subtracted. 

From the practical point of view, the latter algorithm looks simpler. How- 
ever, the fact that UV and IR divergences are not separated causes potential 
problems. In the first approach the expressions 72 — 72 and A( are man- 
ifestly free of IR divergences, so that we might as well set m;, = = 
in their calculation. However, the expression B alone is not, so that in the 
second approach one has to rely on a QCD Ward identity to ensure that 
B is cancelled by other diagrams, and one gets a finite answer. A similar 
statement holds true for the other counterterm diagrams. If C is naively 
calculated, the corresponding diagram (Fig. |3|) is UV and IR divergent, and 
if both divergences are regulated dimensionally, the relevant coefficient is 
left undetermined. Again, an electroweak Ward identity provides the correct 
renormalization of g t , which has been used in ref. |J. 

To conclude, the complete effective-theory expressions are somewhat cum- 
bersome, but they allow a safe diagram-by-diagram analysis. In order to 
simplify the calculations, when additional information is used (symmetry ar- 
guments, explicit expressions for bare parameters), and IR divergences are 
controlled, shortcuts are possible. However, once logarithms have to be re- 
summed, one has to return to the complete expressions developed in the main 
part of the present paper. 
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Figure 1: Sample diagram in a NLO calculation: The coefficient of the 1/e 
pole is denoted by A. 
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Figure 2: Counterterm diagrams: B, C, D, and E are the coefficients of the 
1/e poles. An open square denotes the pole part of the subdiagram; a filled 
square stands for the finite part. 
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Figure 3: Subdiagram for C and E. 
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